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We discuss the CP nature of the Yukawa couplings of the Higgs boson in the framework of a 
complex two Higgs doublet model (C2HDM). After analysing all data gathered during the Large 
Hadron Collider run 1, the measurement of the Higgs couplings to the remaining SM particles already 
restricts the parameter space of many extensions of the SM. However, there is still room for very 
large CP-odd Yukawa couplings to light quarks and leptons while the top-quark Yukawa coupling 
is already very constrained by current data. Although indirect measurements of electric dipole 
moments play a very important role in constraining the pseudo-scalar components of the Yukawa 
couplings, we argue that a direct measurement of the ratio of pseudoscalar to scalar couplings should 
be one of the top priorities for the LHC run 2. 


I. INTRODUCTION 

The Higgs boson discovery by the ATLAS [T] and CMS [ 2 ] collaborations at the Large Hadron Collider (LHC) 
has triggered a number of studies on multi-Higgs extension of the Standard Model (SM). Although the measured 
Higgs couplings show a very good agreement with the SM predictions there is still room for interesting non-SM 
features to be explored at the next LHC run. In fact, many multi-Higgs models provide interesting scenarios 
as is the case of the complex two-Higgs double model (C2HDM). The 2HDM was proposed by T. D. Lee [3] 
as a means to explain the matter-antimatter asymmetry of the universe by allowing for an extra source of 
CP-violation in the potential (see laig for a review). The existing experimental data and in particular the one 
recently analysed at the LHC has been used in several studies with the goal of constraining the parameter space 
of the C 2 HDM [BHS] or just the Yukawa couplings m- 

In this work we analyse C 2 HDM scenarios where the scalar component of the SM-like Higgs Yukawa couplings 
to down-type quarks or to leptons vanish. The corresponding CP-odd component has to be non-zero for the 
model to be in agreement with the LHC results. We will also discuss situations where the Yukawa coupling is 
shared by the CP-even and CP-odd components of the Higgs boson. Our approach is driven both by the current 
measurements and by the predictions for the next LHC run. The processes pp ^ h ^ WW{ZZ), pp ^ h ^ 
and pp ^ h ^ are at present measured with an accuracy of about 20%. The expected accuracies for the 

signal strengths of different Higgs decay modes were presented by the ATLAS m and CMS [ 12 ] collaborations 
(see also [13]) for ^/s = 14 TeV and for 300 and 3000 fb~^ of integrated luminosities. The predictions for the 
signal strengths with the final states VV, 77 and r+r“ will be used to understand how the model will perform 
at the end of the next LHC run because as shown in [9] they reproduce quantitatively the effect of all possible 
final states in the Higgs decay. Therefore, the predicted accuracies for the signal strength lead us to consider 
situations where, at 13 TeV, the rates are measured within either 10% or 5% of the SM prediction. We note 
that there is no visible difference in the plots when the energy is changed from 13 to 14 TeV as discussed in [9]. 
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II. THE C2HDM 


The allowed parameter space of the C 2 HDM was recently reviewed in [14] (see also j6ll9l [T5H2Q] h In this section 
we will briefly describe the C2HDM, a complex 2HDM with a softly broken Z 2 symmetry (fi (f 1,(^2 ~ 02 - 

We write the scalar potential as [5| 

Vh = +ml2\(p2?-m\2(p\(p2-{rn\2)* 

All / |4 , A 2 I 


+ + A 3 I 01 PI 02 P + A 4 (0102) (020l) 

+ ^(0102)^ + ^(020l)^ 5 


( 1 ) 


where all couplings except mi 2 and A 5 are real due to the hermiticity of the potential and arg(A 5 ) 7 ^ 2 arg(mi 2 ), 
so that the two phases cannot be removed simultaneously [T5] . 

We work in a basis where the vacuum expectation values (vevs) are real. The corresponding CP-conserving 
2HDM, is obtained from the C2HDM by taking mi 2 and A 5 real. Defining the scalar doublets as 


Ti 


Pt 


V \-^{V 2 + r ]2 + 1 x 2 ) 

with V = V 2 = (\/2G^)“^/^ = 246 GeV, they can be written in the Higgs basis as EH [ 22 ] 


w 

H 2 


_ ( cp sp 
-Sg Cg 


01 

02 


( 2 ) 


(3) 


where tan0 = V 2 1vi, cg = cos0, and sg = sin0. In the Higgs basis the second doublet does not get a vev and 
the Goldstone bosons are in the first doublet. 

Defining 773 as the neutral imaginary component of the H 2 doublet, the mass eigenstates are obtained from 
the three neutral states via the rotation matrix R 



= R 



(4) 


which will diagonalize the mass matrix of the neutral states via 

RM^ R^ = diag m\) , (5) 

and mi < m 2 < m 3 are the masses of the neutral Higgs particles. We parametrize the mixing matrix R as m 


C 1 C 2 S 1 C 2 S2 

i?=| -(C 1 S 2 S 3 + S 1 C 3 ) C 1 C 3 -S 1 S 253 C 2 S 3 

—C1S2C3 + 5153 —(C1S3 + S1S2C3) C2C3 


with Si = sino^i and q = cosa^ {i = I, 2, 3) and 

— 7r/2 < ai < 7r/2, —7r/2 < a2 < 7r/2, 


—7r/2 < 0^3 < 7r/2. 


(6) 


(7) 


We choose the 9 independent parameters of the C2HDM to be v, tan0, m^±, o^i, 0 ^ 2 , as, mi, m 2 , and 
Re(mi 2 ). The mass of the heavier neutral scalar is a dependent parameter given by 

miRi3(Ri2tan0 - Ru) P ml R23(^22tan0 - R21) 


mo = 


( 8 ) 


RssiRsi — R 32 tan 0 ) 

The parameter space will be constrained by the condition m 3 > m 2 . 

The Higgs coupling to gauge bosons is [ 6 ] 

C = CgRii P SgRi2 = COS ((^ 2 ) COS (o^i — 0) . (9) 

Regarding the Yukawa couplings, the Z 2 symmetry is extended to the Yukawa Lagrangian [23| to avoid flavour 
changing neutral currents (FGNG). The up-type quarks couple to 02 and the usual four models are obtained by 
coupling down-type quarks and charged leptons to 02 (Type I) or to 0i (Type H); or by coupling the down-type 
quarks to 0i and the charged leptons to 02 (Flipped) or finally by coupling the down-type quarks to 02 and the 
charged leptons to 0i (Lepton Specific). The Yukawa couplings can then be written, relative to the SM ones, 
as a + ib^s with the coefficients presented in table [l| 
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Type I 


Type II 

Lepton 

Specihc 

Flipped 

Up 


Ri2 

S/3 

-*C/3^75 



Down 


Rii 




Leptons 


C/3 





TABLE I: Yukawa couplings of the lightest scalar, hi, in the form a + ib'y^- 


III. RESULTS AND DISCUSSION 

In order to perform our analysis we generate points in parameter space in the following intervals: the lightest 
neutral scalar is mi = 125 GeV the angles <^ 1 , 2,3 all vary in the interval [—7r/2, 7r/2], 1 < tan^d < 30, 
mi < m 2 < 900GeV and —(900GeV)^ < Re(mi 2 ) < (900GeV)^. The points are generated randomly subject 
to the following constraints, 

• B-physics - b ^ sj, in Type II/F we choose the range for the charged Higgs mass as 340 GeV < mH± < 
900 GeV [26], while in Type I/LS the range is 100 GeV < mH± < 900 GeV. The remaining constraints 
from B-physics [27l[28] (and from the Rij = T{Z ^ bb)/T{Z hadrons) [29] measurement) force tan^d ^ 1 
for all models; 

• LEP - The charged Higgs mass is above 100 GeV due to LEP searches on e+e“ ^ [30] (we also 

consider the LHG results on pp ^ tt{-^ H^b) [3T||32]). Very light neutral scalars are also constrained by 
LEP results [33] : 

• LHC - bounds on heavy scalars - The most relevant searches for the C2HDM are pp ^ (j) ^ 
W+W-{ZZ) [31135] and pp ^ (j) ^ T+T- IMIEZI, where 0 is a spin zero particle; 

• Theoretical constraints - the potential is bounded from below [38], perturbative unitarity is en¬ 
forced [39ll4T] and all allowed points conform to the oblique radiative parameters [421144] . 

Finally we consider the results stemming from the the 125 GeV Higgs couplings measurements. The signal 
strength is defined as 


hi ^ crBRjhi -)> /) 


( 10 ) 


where a is the Higgs boson production cross section and BR(/i^ ^ /) is the branching ratio of the hi decay 
into the final state /; and BW^{h f) are the respective quantities calculated in the SM. Values for 
the cross sections were obtained from: HIGLU [45] - gluon fusion at NNLO, together with the corresponding 
expressions for the CP-violating model in [9]; SusHi [46] - bb ^ h nt NNLO; m - Vh (associated production), 
tth and VV h (vector boson fusion). As previously discussed, we will force /iyy, and pirr to be within 
20% of the expected SM value, which at present roughly matches the average precision at Icr. Taking all other 
processes into account has no significant impact on the results as shown in [9] . 

In the G2HDM there is only one way to obtain pure scalar states. When we set 52 = 0 we get R 13 = 0 and 
all pseudoscalar components vanish. However, depending on the model type, there are in principle two ways 
to obtain a vanishing scalar component. One is by setting R 12 = 0. However, as shown in figure (middle), 
values of R 12 ~ 0 are excluded when all constraints are taken into account. The other possibility is to have 
Rii = Cl C2 = 0 which is still allowed as shown in figure (left). This can be obtained by setting either C2 = 0 
or Cl = 0. C 2 = 0 is excluded as it would mean gn^vv = O 5 where V is a massive vector boson. Finally we can 
choose Cl = 0. The values for ap and bp {F = U, D, L) in this scenario are presented in Table [H] 


^ The latest results on the measurement of the Higgs mass are 125.36 ± 0.37 GeV from ATLAS |24| and 125.02 + 0.26 — 0.27 (stat) 
+0.14 — 0.15 (syst) GeV from GMS [25]. 









Toyama International Workshop on Higgs as a Probe of New Physies 2015, 11-15, February, 2015 


4 



FIG. 1: Top: tan/3 as a function of Rn (left), R 12 (middle) and Ri^ for Type 1. Bottom: same but for Type IL The 
rates are taken to be within 20% of the SM predictions. The colours are superimposed with cyan/light-grey for jivv^ 
blue/black for yrr and finally red/dark-grey for with ^/s = 8 TeV. 


Type I au — an — ul — ^ bu — —bo = —bh = — 

Type II an = ul = bn = bL = -S 2 tp 

Type F an = 0 bo = —S 2 tp 

Type LS ul = 0 bh — —S 2 tp 


TABLE II: ar and bp limits for ci = 0 for the four model types. 

As shown in table |n| the scenarios where the scalar component vanishes arise only in models Type II, F and 
LS. In Type II one can have = a l = 0 while in F (LS) only ai:) = 0 (ai, = 0) is possible. In this scenario the 
coupling to gauge bosons is 

C^=slcl. ( 11 ) 

Note that even if 52 = 0 the pseudoscalar component can still be large due to a large value of tan (I. 

We will now discuss in detail the allowed parameter space in the {aF^bp) plane for the different model types. 
We will plot sgn(C) ap (sgn(C) bp) instead of ap (bp) with F = U,D, L, to avoid the dependence on the phase 
conventions in choosing the range for the angles ai. In the left panel of figure we show bp = bp diS a function 
of ap = ap for Type II and y/s = 13 TeV with all rates at 10% (blue/black), 5% (red/dark-grey), and 1% 
(cyan/light-grey). We start by noting that this scenario is still possible with the rates within 5% of the SM 
value at the LHC and at =13 TeV. This scenario can only be excluded by a measurement of the rates if 
the accuracy reaches about 1%. The constraints on the model force \bp\ 1 when \ap\ 0. When \bp\ ~ 1, 
the couplings of the up-type quarks to the lightest Higgs have the form 

4 = (1 - 4) = (1 - n4), bi = 4 = 1A4 (12) 

while the coupling to massive gauge bosons is now 

= {t} - + 1) = (1 - sl)/{l + 4 ) . (13) 

In the right panel of figurewe show bp as a function of ap for Type H with the same colour code. We conclude 
from the plot that the constraint on the values of {ap, bp) are already quite strong and will be much stronger 
in the future just taking into account the measurement of the rates. 
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FIG. 2: Left: sgn(C) bo — sgn(C) bL as a function of sgn(C) ud = sgn(C') ul for Type II and a center of mass energy of 
13 TeV with all rates at 10% (blue/black), 5% (red/dark-grey), and 1% (cyan/light-grey). Right: same, but for sgn(C) 
bu as a function of sgn(C) au- 
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FIG. 3: Left: sgn(C) bu as a function of sgn(C') au for Type I and a center of mass energy of 13 TeV with all rates at 
10% (blue/black) and 5% (red/dark-grey). Right: sgn(C) bu as a function of sgn(C) au for LS and a center of mass 
energy of 13 TeV with all rates at 10% (blue/black) and 5% (red/dark-grey). 


In the left panel of figure]^ we show bu as a function of au for Type I and ^/s = 13 TeV with all rates at 
10% (blue/black) and 5% (red/dark-grey). We should point out that even at 10% there are still allowed points 
close to (a, 6) = (0.5, 0.6) with no dramatic changes occurring for an increase in accuracy to 5%. In the right 
panel we present 6^ as a function of ai, for Type LS with the same colour code. Here again the (au^bu) = (0,1) 
scenario is still allowed with both 10% and 5% accuracy. However, as was previously shown, the wrong sign 
limit is not allowed for the LS model [48l[49]. Nevertheless, in the C2HDM, the scalar component sgii{C) au 
can reach values close to —0.8. Finally, for the up-type and down-type quarks, the plots are very similar to the 
one in the right panel of figure for Type IT 

In the left panel of figure we present the allowed space in the sin(a2-tan/3 plane, for Type H and ^/s = 13 
TeV. Rates at 10% are shown in blue/black while in red/dark-grey we present the points with \ai)\ < 0.1 and 
II^dI — 1| < 0.1 and in green Ibul < 0.05 and Waul — 1| < 0.05. The right panel now shows the allowed space 
in the sin (a 2 -cos ai plane. The purpose of these plots is to pinpoint the main differences between the SM-like 
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FIG. 4: Left: tan^^ as a function of sin 0^2 for Type II and a center of mass energy of 13 TeV with all rates at 10% 
(blue/black). In red/dark-grey we show the points with \aD\ < 0.1 and \\bD\ — 1| < 0.1 and in green |5 d| < 0.05 and 
\\aD \ — 1| < 0.05 Right: same, with tan/3 replaced by cosai. 


scenario, where {\aD\, {bnl) ~ (I 5 O) and the pseudoscalar scenario where {\aD\^ \^d\) ~ la the SM-like 

scenario sin 0^2 ~ 0, tan/3 is not constrained and the allowed values of sin 0^2 grow with increasing coso^i. In the 
pseudoscalar scenario coso^i ~ 0, sin 0^2 and tan^^ are strongly correlated and tan/3 has to be above ~ 3. We 
note that all values of an and bo are allowed provided P bjj ^ 1. 


A. Direct measurements of the CP-violating angle 

We have seen that the precise measurements of the Higgs couplings allows us to constrain both the scalar 
and the pseudoscalar Yukawas in the C2HDM. However, a direct measurement of the relative size of the 
pseudoscalar to scalar coupling is important because it directly probes the Higgs couplings to light quarks and 
leptons. Moreover, when combined with EDMs it can provide universality tests for the CP-odd components of 
the Yukawas. 

The angle that measures the pseudoscalar to scalar ratio, 0^, is defined by 

tan0^ = bilai i = U, D, L , (14) 

and could in principle be measured for all Yukawa couplings. Direct measurements of this ratio in the up-quark 
sector, bu jau, was first proposed in [50] and more recently in [5TH53] . The process pp hjj [54| also allows 
to probe the same vertex as discussed in (55] |56]. In reference [56] an exclusion of > 40° {(ft > 25°) for a 
luminosity of 50 fb“^ (300 fb“^) was obtained for 14 TeV and assuming = 0 as the null hypothesis. A study 
of the T^T~h vertex was proposed in m (see also [SlISS]) and a detailed study taking into account the main 
backgrounds m lead to an estimate in the precision of of 27° (14.3°) for a luminosity of 150 fb ^ (500 
fb“^) and ^/s = 14 Tev. 

The number of independent measurements of ft one needs depends on the C2HDM Yukawa type. For Type 
I one process is enough since <fu = fo = fr- For all other types we need two independent measurements. For 
type H and LS the planned measurements of (ft and fr would be enough while for type F we would need fi,. 
Incompatibility in the measured values of ft and f^ would exclude both Type I and F. 

We will now discuss the behaviour of coso^i (figure left) and sgn{C) ao (figure right) as a function of 
fo = fr for Type IT In figure]^ = 13 TeV and all rates are taken at 10% (blue/black); in red/dark-grey we 
show the points with \aD\ < 0.1 and ||6i:)| —1| <0.1 and in green {bol < 0.05 and ||ai:)| —1| < 0.05. The SM-like 
scenario sgn((7) {ao, bjj) = (1,0) is easily distinguishable from the (0,1) scenario. In fact, a measurement of 
fr even if not very precise would easily exclude one of the scenarios. Obviously, all other scenarios in between 
these two will need more precision (and other measurements) to find the values of scalar and pseudoscalar 
components. The angle is related to a 2 as 


iceafr = —sgjc\ tan 0^2 


tan 0^2 = —cifsg tanfr 


(15) 
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The zero scalar 


limit Type II SM-like limit t „|| 



0T = tan-‘'(bo/ac») (defl) 0^ - 0 ^ tan*‘'(bo/ao) (deq) 

FIG. 5: Left: cosai as a function of ta,n~^{bD/aD) for Type II and a center of mass energy of 13 TeV with all rates at 
10% (blue/black). In red/dark-grey we show the points with \aD\ < 0.1 and \\bD\ — 1| < 0.1 and in green \bD\ < 0.05 
and Waol — 1| < 0.05. Right: same, with cosai replaced by sgn(C)ai:). 

and therefore a measurement of the angle (j)r does not directly constrain the angle a 2 but rather a relation 
between the three angles. A measurement of (ft and would give us two independent relations to determine 
the three angles. 


B. Constraints from EDM 

The C2HDM, as all models with CP violation, are constrained by bounds arising from the measurement of 
electric dipole moments (EDMs) of neutrons, atoms and molecules. The parameter space of the C2HDM was 
analysed in lailolISlHel] and ref. [7] found that the most stringent bounds are obtained using the results from 
the ACME Collaboration [65], except when cancellations among the neutral scalars occur. These cancellations 
were pointed out in [63[ [64] and arise due to orthogonality of the R matrix in the case of almost degenerate 
scalars |9]. 

The scenarios we discuss have the couplings of the up-type sector (top quark) very close to the SM ones. 
Indeed, only the couplings in the down-type sector, namely tan lepton and the b-quark Yukawa couplings, 
are still allowed by data to have a vanishing scalar component. Due to the universality of the lepton Yukawa 
couplings the electron EDM also restricts the tan Yukawa. In a type II model this in turn also restricts the 
b-quark Yukawa of the SM-like Higgs. However, this is completely irrelevant in the Elipped model, where the 
charged leptons couple as the up-type quarks. We have shown in El that in a preliminary scan over the 
parameter space we have found points which pass all constraints including ACME’s. 

A dedicated study of the EDM contributions in the Type H C2HDM, where there are several sources of CP 
violation and where the partial cancellations of the various scalars is dully taken into account is in progress [66] . 
In addition, one should keep in mind that, as pointed out in ref. ElEZl, the future bounds from the EDMs 
can have a strong impact on the C2HDM. In the future, the interplay between the EDM bounds and the data 
from the LHC Run 2 will pose relevant new constraints in the complex 2HDM in general, and in particular for 
the scenarios presented in this work. 


IV. CONCLUSIONS 

We have discussed the interesting possibility of having a vanishing scalar component is some of the Yukawa 
couplings, namely the couplings of the lightest Higgs to down-type quarks and/or to leptons. These scenarios 
can occur for Type H, E and LS, and the pseudoscalar component plays the role of the scalar component in 
assuring the measured rates at the LHC. A direct measurement of the angles that gauge the ratio of pseudoscalar 
to scalar components is needed to further constrain the model. In particular, the measurement of 0^-, the angle 
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for the T^T~h vertex, will allow to either confirm or to rule out the scenario of a vanishing scalar, even with a 
poor accuracy. We have also noted that for the Type F, only a direct measurement oi (j)D in a process involving 
the bbh vertex would probe the vanishing scalar scenario. Finally a future linear collider [68] |69] will certainly 
help to further probe the vanishing scalar scenarios. 
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